Temperature Dependence of Upper Critical Field as an Indicator of boson Effects in 

Superconductivity in Nd2_ x Ce a Cu04_ 



-v 



V. F. Gantmakhcr *, G. A. Emel'chcnko, I. G. Naumcnko, and G. E. Tsydynzhapov 
Institute of Solid State Physics RAS, 142432 Chernogolovka, Russia 



o 
o 

(N 

c 

Hi 



G ■ 
O ■ 

■ 

1 1 
5h ■ 
Oh 

3 : 

. 



I 

C 

o 
o 



> 

00 
(N 

o 



i 

c 

o 
o 



The temperature dependence of upper critical field B C 2 
was determined from the shift of resistive transition AT(B) 
in nearly optimally doped Nd2- a; Ce 2 :Cu04_j / single crystals. 
Within the experimental accuracy, the weak-field data are 
described by power function B C 2 oc (AT) 3 / 2 . This result is 
compared with the data on heat capacity and analyzed in the 
context of possible manifestations of boson effects in super- 
conductivity. The T dependence of B C 2 persists down to the 
lowest temperatures, but the numerical values of B C 2 below 1 
K are different for different samples. 

PACS numbers: 74.20.Mn, 74.25.Dw, 74.72.-h.22 

There are grounds to believe that high-temperature su- 
perconductivity (HTSC) is not described by the BCS the- 
ory. One of them consists in the relationship between the 
density n of Cooper pairs and the coherence length £ (the 
pair size). In HTSC cuprates, superconductivity is due 
to the carriers in the Cu02 plane. Like in all 2D systems, 
the density of states gp at the Fermi level in the Cu02 
plane does not depend on the carrier concentration in the 
normal state and, according to measurements, is equal to 
gp = 2.5 ■ 10 K _1 per one structural unit of Cu02 (this 
value is nearly the same for all cuprate families, see, e.g., 
[1], Ch. 13). Assuming that the superconducting gap A 
is of the order of transition temperature T c , one estimates 
the mean distance r = rT 1 ! 2 ps (gpA)~ 1 / 2 between the 
pairs in the Cu0 2 plane at 25 A for T c ps 100 K and 75 A 
for T c psIOK. These r values should be compared with 
the typical coherence length £ ps 20 A in the ab plane 
[1], so that r > £ in the HTSC materials. Inasmuch as 
the BCS theory introduces Cooper pairs to describe the 
Fermi-liquid ground state as a whole, its validity for the 
description of HTSC is not obvious. This causes inter- 
est in the models of Bose superconductivity for which 
r 3> £ an d which are based on Bose-Einstein conden- 
sation (BEC) in a system of charged bosons [2-4]. The 
experimental evidences for the boson effects in HTSC are 
presently intensively accumulated. 

One such evidence can be expected to obtain from the 
measurements of the temperature dependence of mag- 
netic field B C 2 destroying superconductivity. In the BCS 
theory, the B C 2(T/T C )/B C 2(0) function is linear in the 
vicinity of T/T c = 1; it monotonically increases to sat- 
uration near the zero temperature and almost coincides 
with the limiting value even at T/T c = 0.2 [5]. However, 
in most cases, the HTSC materials behave in a differ- 
ent manner and demonstrate positive second derivative 
d 2 B C 2/dT 2 over the entire temperature range. 

The B C 2{T) measurements are mainly based on an 



analysis of the resistive transition. Two types of be- 
havior are known for the resistive transition of HTSC 
materials in a magnetic field. For one of them, the tran- 
sition is sizably broadened in a magnetic field, so that 
it is hard and even practically impossible to gain from 
it any information about the B&iT) dependence. The 
other transition is shifted in a magnetic field to lower 
temperatures and either remains undistorted, as in usual 
superconductors, or undergoes an insignificant distortion. 
This usually occurs for those members of HTSC families 
in which T c < 20 K. The transition shift in these mate- 
rials is naturally explained by the field-induced destruc- 
tion of superconductivity. Irrespective of the mechanism 
of dissipative processes in the superconducting state, the 
spectrum rearrangement and the appearance of super- 
conducting pairing should necessarily affect the R(T) re- 
sistance. With this preposition, one can readily construct 
the B c2 (T) function. 

Almost in all HTSC cuprates such as the Tl-based 
[6] and Bi- based [7] families and the LaSrCuO [8] and 
Nd(Sm)CeCuO [9-11] families, as well as in the Zn-doped 
[12] or oxygen-deficient [13] YBaCuO, the B c2 {T) func- 
tion derived from the shift of resistive transition has the 
positive second derivative over the whole temperature 
range < T/T c < 1 and shows a tendency to diverge 
at small T/T c values. Most discussion over the B C 2(T) 
curves concentrated precisely on this divergence and con- 
sidered it as the most dramatic departure from the BCS 
theory. At the same time, the behavior of the B C 2(T) 
function near T c is also quite informative. Contrary to 
expectations, almost in all cases where the field-induced 
resistive-transition shift in HTSC cuprates pro-ceeds in a 
parallel manner, the experimental data indicate that the 
0B C 2/dT derivative is zero at the T c point [6-13]. 

The dB c /dT derivative of critical field in the T c point 
is related to the free energy F and heat capacity C in 
this point by the well-known Rutgers formula: 
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Inasmuch as the thermodynamic critical field B c is differ- 
ent from the upper critical field B&, Eq. (1) can be used 
only for qualitative estimates. However, being based on 
thermodynamics, this equation is very useful. 

In usual superconductors, F s — F n cx (T c — T) 2 , so that 
the heat capacity undergoes a jump and B c is linear in 
(T c - T). In the BEC case, F s - F n cx (T c - T) 3 , so 
that the heat capacity is a continuous function in the 
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transition point [14]. It then immediately follows that 
dB c /dT = and 

B c oc (T c - Tf' 2 . (2) 

Of course, one can hardly imagine that the Fermi gas 
suddenly and completely transforms into a Bose gas at 
low temperatures. It was assumed in [4] that bosons 
appear in small pockets of the k-space near the Fermi 
level. In the isotropic model, one can only speak about 
pairing of sufficiently energetic fermions, as in the BCS 
theory. This kind of model has been proposed in [15]. 
Nevertheless, Eq. (2) deserves a serious experimental 
verification. 

Such was the motivation of our work consisting in the 
measurement and analysis of the field-induced shift of 
resistive transition in Nd2-xCe a; Cu04_ a single crystals. 
We will discuss separately the behavior of the B c2 field 
in the vicinity of T c and at low temperatures. 

Experiment. (NdCe)2 CUO4 single crystals were 
grown from a mixture of components taken in the mo- 
lar ratio Nd2 3 :Ce0 2 :CuO = 1:0.05:11 in a crucible 
made from yttrium-stabilized zirconium dioxide. The 
use of a modified growth regime markedly reduced the 
time of interaction between the melt and the crucible at 
high temperature. Owing to the accelerated-decelerated 
crucible rotation, the melt was intensively stirred so that 
the homogcnization time for the molten solution did not 
exceed 1 h at a temperature near 1150 °C. The growth 
was carried out for several hours upon slow cooling (6 
K/h) under the conditions of morphologically stable crys- 
tallization front (dT/dx > 10 K/cm), after which the 
crucible was decanted and cooled at a rate of 30-50 
K/h to ambient temperature. The crystals were shaped 
like platelets of thickness 20-40 mm. Their composition 
Ndi^Ceo.isCuOz was determined by local X-ray spec- 
troscopic analysis. The analysis revealed Zn traces in the 
crystals at a level of 0.1 wt %. Initially the crystals did 
not show superconducting transition above 4.2 K. The 
superconducting transition at T c w 20 K appeared after 
15-h annealing at 900 °C in an argon atmosphere. 

Measurements were made for two plates approximately 
1x2 mm in size. The silver paste contacts were fused 
in the air at a temperature of 350 °C. Four contacts in 
sample 1 were arranged 0.5 mm apart in a row on one 
side of the plate. The potential contacts in sample 2 
were placed on the opposite side of the plate beneath the 
current contacts, allowing the measuring current to be 
directed both along and transversely the ab plane. This 
did not affect the results. The resistance was measured 
by the standard method using a lock-in nanovoltmeter at 
a frequency of 13 Hz. The measuring current was small 
enough for the linear regime and the absence of overheat- 
ing to be provided down to the lowest temperatures. The 
magnetic field was directed along the normal to the plate 
(c axis) . Measurements were performed over the temper- 
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FIG. 1. The R(T) curves for sample 2 in magnetic fields 
(from right to left) from to 7 kOe. The dashed lines are the 
straight line R n (T) and the straight lines at the levels of 0.67, 
0.5, and 0.2 of R n (T). The method of determining the onset 
of transition is demonstrated and the T c ; fields from which 
the shifts were measured are shown 

ature range from 25 K to 25 mK. [] The onset of zero-field 
superconducting transition in both samples occurred at 
about 20.5 K. 

The measurements gave identical results for both sam- 
ples. Figure 1 demonstrates a scries of low-field R(T) 
curves for sample 2. At high temperatures, all curves 
show the same asymptotic behavior R n (T) above the 
transition, and one can assume that the R n function 
does not depend on B at T > 10-12 K. The zero-field 
transition shows a certain structure, which, however, is 
smoothed out even at 100-200 Oe. The field effect mainly 
amounts to shifting the transition to lower temperatures. 
The degree in which this shift is parallel can be checked 
by comparing the shift of the onset of transition with 
the shifts of the R(T) function at different levels: 0.2R n , 
0.5i?„, and 0.67i?„ (see curves in Fig. 1). For the paral- 
lel shift, all constructions in Fig. 1 should give the same 
function B c2 (AT) 1 where AT = T ci - T and T cl is the 
temperature corresponding to the same level on the ini- 
tial curve R{T,B = 0). The log-log plots of the shifts 
are shown by different symbols in Fig. 2a for all four lev- 
els. The systematic deviations of the symbols from the 



x The low-temperature measurements in strong magnetic 
fields were carried out at the NHMFL (Tallahassee, Fla., 
USA). 
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FIG. 2. (a) Plots of the field vs. shift at different levels in 
this field; (b) the same for the averaged shifts for two samples. 

straight line 

B c2 = (AT)' 3 , (3) 

constructed by averaging the results for all points are 
small for each of the symbols. This implies that the dis- 
tortions of the transition shape are small as compared 
to its shift. The scatter of points in low fields is mainly 
caused by the fine structure of the R(T, B = 0) curve, 
which serves as a reference in the determination of the 
shift AT. 

The coefficient (3 was determined from the slope of the 
straight line passing through the averaged AT shifts (Fig. 
2b). The curve processing for sample 2 (Fig. 1) yields 
/? w 1.4, and the processing of analogous curves for sam- 
ple 1 yields (3 « 1.5. 

The resistances for both crystals decreased in a rela- 
tively narrow temperature range not to zero; one can see 
in Fig. 1 that, starting at the level of ~ 0.1, a slanting 
tail appears. The same tail for sample 1 starts at a higher 
level of ~ 0.2. In this work, we will analyze only the up- 
per portion of the transition, assuming that the electron 
spectrum is rearranged into the form characteristic of the 
superconducting state precisely in this region. 

Figure 3 shows the R(B) functions for very low tem- 
peratures T/T c0 < 0.05. In this region, the normal resis- 
tance depends, though weakly, on a magnetic field, while 
the onset of transition is clearly defined and its shift is 
easily detected even upon changing temperature below 
T/T c0 = 0.005. When considering the B c2 {T) functions 
in this region (see inset in Fig. 3), two fact are notewor- 
thy. First, B c2 does not show tendency to diverge near 
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FIG. 3. The R(B) curves for sample 1 at temperatures 
(from left to right) from 0.5 K to 25 mK. Inset: the field of 
the onset of transition at low temperatures for both samples. 

zero temperature; although the derivative of B C 2{T) is 
large below 0.5 K, the function is linear within the ex- 
perimental accuracy and extrapolated to a finite value 
-B C 2(0) (similar result was obtained previously for thal- 
lium crystals [6]). Second, the critical fields at low tem- 
perature are equal to 69 and 80 kOe for samples 1 and 
2, respectively, i.e., differ by more than 10%, inspite of 
the fact that the crystals were from the same batch and 
their T c values coincide. 

The graph of B c2 (T) over the entire temperature range 
is shown in the inset in Fig. 4; as in other HTSC cuprates, 
the second derivative d 2 B c2 /dT 2 > for all temperatures 
(cf., e.g., [6, 7]). 

Discussion. It follows from the preceding section that 
our data for the vicinity of T c are consistent, within the 
experimental accuracy, with Eq. (2). It would have been 
instructive to compare these data with the data on heat 
capacity, but, unfortunately, in the works where the heat 
capacity of Nd2- x Ce a: Cu04_ y was measured [16] the con- 
tribution of critical fluctuations near T c was not deter- 
mined. Nevertheless, it is known that the measurements 
of heat capacity of the HTSC materials show strong dis- 
similarity from usual superconductors [17] but do not al- 
low the discrimination between the BCS and BEC mod- 
els. These problems can be illustrated by comparing the 
results of measurements of the resistance and heat capac- 
ity of the thallium high-T c superconductor. No explicit 
jump in heat capacity is observed for this compound even 
at zero field, although the contribution from the criti- 
cal fluctuations is undoubtedly present in the tempera- 
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FIG. 4. Comparison of the experimental B C 2 values for 
sample 1 with Eqs. (2) and (4). Inset: the B C 2(T) function 
for sample 2 in the whole temperature region. The line is a 
guide to the eye. 

ture range 16-10 K [18]; this contribution is reduced by 
approximately one-half in a field of 0.4 T and remains 
virtually in the same temperature range. At the same 
time, the resistive measurements by the same experimen- 
tal group [6] suggest that a field of 0.4 T shifts the tran- 
sition by 25% from 16 to 12 K. 

In connection with this contradiction, an interesting re- 
mark was made in [19], where the numerical calculations 
were carried out for the heat capacity of an ideal charged 
Bose gas in a weak magnetic field. It is well-known that 
the BEC does not occur in the ideal charged Bose gas 
in a uniform magnetic field [20] because the density of 
states diverges at the lower Landau level of the spectra 
of charged bosons. This implies that the transition oc- 
curs only at an isolated point in the (T, B) plane. The 
magnetic field in this plane is scaled by the comparison of 
the cyclotron energy TieB/mc with T c . Substituting the 
free electron charge and mass for e and to, respectively, 
one arrives at the value of 8 T for the characteristic field 
at T c — 16 K. On this scale, the above-mentioned field 
of 0.4 T is as small as 0.05. As long as the field is low, 
the phase trajectory again passes through the vicinity of 
the transition point in the (T, B) plane upon changing 
T, but, as the field increases, the "impact parameter" in- 
creases, while the contribution of critical fluctuations de- 
creases. However, the temperature interval correspond- 
ing to the small impact parameters does not change. In 
the case that the transition is the BEC in a weakly non- 



ideal charged Bose gas, this contribution is hidden from 
view at the lower temperature where the transition oc- 
curs in magnetic field. Then, strange as it may seem, the 
resistive measurements provide the more reliable infor- 
mation on the transition position than the heat capacity 
measurements do. 

According to the results obtained for the immediate 
vicinity of T c , the behavior of the B C 2 (T) function should 
be compared with the predictions of the superconductiv- 
ity models in a nonidcal Bose gas. Due to the boson scat- 
tering by impurities or to the boson-boson interaction, 
the critical field in a weakly nonideal Bose gas behaves 
as [21] 

B c2 oct- a (l-i 3 / 2 ) 3 / 2 , t = T/T c , (4) 

where, depending on the particular model, the exponent 
a is equal to 1 or 3/2 [21, 22]. At t ->■ 1, function (4) 
takes the asymptotic form (2). It is seen in Fig. 4 that 
the experimental points deviate in the proper direction 
from the asymptote and, on the whole, correspond well to 
Eq. (4). A more detailed comparison is hardly pertinent, 
as long as the theories [21, 22] do not allow for the field- 
induced pair decay into fermions. 

Conclusions. The field-induced distortion of the 
shape of resistive superconducting transition in the 
Nd2- rr Cc a ;Cu04_j ) single crystals is appreciably smaller 
than the transition shift. This allows the measurement of 
the B c2 (T) function. As zero-field T c is approached, the 
B C 2 field behaves as a power function B C 2 oc (AT)@ with 
f3 w 1.5 and, correspondingly, with a horizontal tangent 
dB c2 /dT = 0. This should imply the absence of jump in 
heat capacity at the zero-field phase transition. Such a 
behavior is precisely that which is expected for the heat 
capacity and critical field in the BEC of a charged Bose 
gas. For this reason, one of the possible conclusions that 
can be drawn from such a behavior of B c2 (T) near T c is 
that the description of superconductivity of the HTSC 
materials should involve the BEC elements, i.e., make al- 
lowance for the fact that fermions near the Fermi level 
tend to form bosons at temperatures above T c . The T 
dependence of B c2 persists down to the lowest tempera- 
tures, although, probably, the B c2 values in this region 
depend on lattice defects. 
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